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ABSTRACT 

Ram pressure stripping, i.e. the removal of a galaxy's gas disk due to its motion through the intracluster medium of a galaxy cluster, 
appears to be a common phenomenon. Not every galaxy, however, is completely stripped of its gas disk. If the ram pressure is 
insufficiently strong, only the outer parts of the gas disk are removed, and t he inner gas disk is retained by the galaxy. One example 
of such a case is the Virgo spiral NGC 4402. Observations of NGC 4402 dCrowl et alj|2005t) reveal structures at the leading edge 
of the gas disk, which resemble the characteristic finger-like structures produced by the Rayleigh-Taylor (RT) instability. We argue, 
however, that the RT instability is unlikely to be responsible for these structures. We demonstrate that the conditions under which a 
galaxy's disk gas experiences ram pressure stripping are identical to those that lead to RT instability. If the galaxy's gravity prevents 
ram pressure stripping of the inner disk, it also prevents the RT instability. In contrast, the stripped gas could still be subject to RT 
instability, and we discuss consequences for the stripped gas. 
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1. Introduction 

Differences between field and cluster galaxies imply that the 
evolution of galaxies is influenced strongly by their environ- 
ments. In addition to the gravitational interaction between the 
cluster galaxies themselves, ram pressure stripping (RPS), the 
(partial) removal of a galaxy's gas due to its motion through the 
intracluster medium (ICM), is thought to play an import ant role 
for the evolution of galaxies (see e.g. Ivan Gorkom 2004). 

RPS has been mod el led nu merically by d i fferen t 
groups (e. g. lAbadi et all |1999t ISchulz & Struckl I2001 



Ouili s et all l2000t IVollmer et all 200 It 



Mar colini et alj 
hereafter RH05; 
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2003t IR oediger & Henslerl I2005I 
Roediger & Briigge n| I2006I here after RB06; IRoediger et al.l 
2006; Roediger & Bruggen[ |2007l hereafter RB07). In environ- 
ments in which ram pressure stripping is prevalent, such as 
cluster centres, the stripping is very effective, but the central 
part of the gas disk can be retained by the galaxy for quite some 
time. 

In principle, the remaining gas disk can be affected by the 
Rayleigh-Taylor (RT) instability, which occurs when a lighter 
fluid supports a heavier one against a gravitational field or an ac- 
celeration. An equivalent situation occurs when a denser fluid 
is accelerated by a lighter one, e.g. when a dense gas cloud 
is exposed to a wind of rarefied gas (or when a gas cloud 
moves through this ambient gas). Such a gas cloud is expected 
to be destroyed by the RT instability, and subsequent Kelvin- 
Helmholtz (KH) instability, unless stabilising mechanisms are 
in place such as (self- )gravity, surface tensio n, magnetic fields 
( Chandrasekharl [l96lh . or heat conduction dVieser & HensleJ 
120071) . 



Thus, as a (disk) galaxy moves through the ICM and its gas is 
stripped by RPS, it should be affected by the RT instability. The 
detection of a filamentary leading surface in NGC 4402, a Virgo 
cluster galaxy with clear signs of RPS, lends support to the idea 
that the R T instability does aff ect the remaining gas disk. Recent 
images bv lCrowl et al.l ((2005) reveal finger-like structures at the 
galaxy's upstream side, which are slightly inclined with respect 
to the rotation axis of the disk. If these pillars are interpreted as 
the results of the RT instability, the direction of motion of the 
galaxy that the y imply would b e identical to that implied by HI 
observations dCrowl et al.ll2.005h . 

Such structures are not, however, seen in the simulations. In 
Sect. 13.11 we demonstrate analytically that the RT instability is 
suppressed in the remaining gas disk due to the gravitation of the 
galaxy. In addition, we discuss the role of the RT instability in 
deciding the fate of stripped gas in Sect. 13.21 To start with, we 
summarise basic analytical estimates concerning the RT insta- 
bility in Sect. [2] 



2. The Rayleigh-Taylor instability 

A detailed derivation of stability condit i ons and instability 
growth rates is provided by Chandrasekhar dl96ll) . 

We consider the simplest case of two superimposed, incom- 
pressible, inviscid, unmagnetised fluids. These fluids are sepa- 
rated by a contact discontinuity that is perpendicular to the ef- 
fective homogeneous (gravitational) acceleration, a. Referring to 
the direction of a as "downwards", this setup is unstable if the 
upper fluid has a higher density than the lower one. Any per- 
turbation at the interface will be amplified by the RT instabil- 
ity, developing characteristic finger-like structures. In the linear 
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Fig. 1. Sketch of a cloud moving through an ambient medium, in 
the rest frame of the ambient medium (left), and in the cloud's 
rest frame (right). The direction of the gravitational acceleration 
at the cloud's upstream side, due to its self-gravity, is shown by 
the vector g. 



regime, small perturbations grow at the most rapid rate. Also a 
continuous stratification will be unstable if the density increases 
upwards. 

Here we are interested in the situation sketched in Fig. Q] 
a dense cloud of mass, M c \, density, p c i, and cross section A is 
moving with velocity, -v am b, through an ambient gas of density 
pamb, or, equivalently, is exposed to a wind of ambient gas with 
velocity v am t,. The relevant acceleration is the one exerted by the 
drag force Fq'- 



a D 



Fp _ 2 



c w Ap amb vl mh 



cl 
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(1) 



where cw is the drag coefficient, and the mean column density 
of the cloud is given by 



£h = 



Mel 

A 



4/3jr/£p c i 



7lR 2 , 
cl 



= 4/3fl cl p cl . 



(2) 



The RT instability will be suppressed, if a e g > 0, or 

8 > a D . (4) 

Another mechanism to stabilise the cloud against RT insta- 
bility, is provided by magnetic fields perpendicular to the cloud 
surface. Magnetic fields act like a surface tension and suppress 
the growth of perturbations of scale-length A < A c with 

B 2 cos 2 6» B 2 cos 2 6>E cl 



Ar = 



y"o(Pcl - Pamb)flD 



(5) 



where 6 is the angle between the magnetic field dire ction and 
the direction of the perturbation ( Chandrasekharl [l96 ll) . We as- 
sume that p c i » pa m b- If A c is larger than the size of the cloud, it 
should therefore be resistant to RT instability because the largest 
possible perturbations have a size that is similar to that of the 
cloud. 

3. Application 

3.1. The remaining gas disk 

We now investigate to what extent the RT instability can develop 
in a disk galaxy that is either moving face-on through the ICM, 
or is exposed to an ICM wind of density picm and velocity vicm- 
In Fig.Q] we replace the cloud by a disk galaxy with gas surface 
density, X gas . The gravitational acceleration due to the galaxy po- 
tential, g, at the upstream side of the galaxy, points downstream. 
According to Eq. |4] the condition for RT instability to occur is 



1 PICM V ICM 



2* 



Picm v icm 



(6) 



(7) 



For a disk galaxy, both, g and £ gas depend on the radial posi- 
tion within the disk and both typically decrease with disk radius, 
r. In the outer parts (where g£ gas is smaller), the RT instability 
can occur, while it is suppressed in the inner part, where g£ gas 
is larger. An identical criterion is applied to decide the radius 
out to which gas will be stripped from the disk by ram-pressure 
(see e.g. RH05, RB06). The ram pressure, jC w picM v icm' can ^ e 
compared to the gravitational restoring force, g£ gas , as written in 
Eq. Q At radii where this inequality applies, the gas will be ram 
pressure stripped; if the inequality does not apply, the expression 
on its left-hand side will dominate, and the galaxy retains its gas. 
In this second case, the inner gas disk is not stripped and is pro- 
tected against RT instability by the gravity of the galaxy. This 
behaviour is reproduced by simulations (see e.g. Fig. |2j. 

We note that tangential magnetic fields would stabilise the 
remaining disk even further. 



The drag coefficient depends on the cloud's shape and surface 3,2. The stripped gas 



properties, where ^cw ~ 1 may be used as an approximation. 
At the upstream side of the cloud, in the rest frame of the cloud 
the drag force appears as a pseudo force pointing upstream, the 
correct direction to produce the RT instability at the cloud's up- 
stream edge. 

A possible way to stabilise the cloud's surface against the 
RT instability would be a gravitational field associated with the 
cloud. The gravitational acceleration, g, at the cloud's upstream 
edge works in the opposite direction of ap ("inwards" for the 
cloud), so that the relevant acceleration for the development of 
the RT instability is the effective acceleration 



-a D + g 



(3) 



According to the results of the previous section, the galaxy's 
gravitational field cannot prevent the RT instability of the 
stripped material. In the absence of other stabilising effects, the 
stripped gas should therefore fragmentate. The simulations of 
e.g. RH05, RB06 and RB07 neglected all possible stabilising ef- 
fects, and therefore led to fragmentation. Figure [2] displays two 
simulation snapshots that differ only in spatial resolution. The 
resolution fixes the lengthscale of the smallest possible pertur- 
bations. According to analytical estimates, the smallest perturba- 
tions grow most rapidly, thus the stripped gas should fragmentate 
into smaller clouds, in simulations of higher resolution. This be- 
haviour is in fact observed in, both, our 2D and 3D simulations. 
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Fig. 2. Simulation snapshots from simulations described in 
RB07. Identical conditions except for the resolution: in rhs 
panel, the spatial resolution is 2 times better than in the lhs panel. 
The remaining gas disk is not fragmented, but the stripped gas 
is. The fragmentation is stronger for higher resolution. 



field, but can be averted by tangential magnetic fields. We note 
that the relevant surface for KH instability are the sides of the 
gas clouds. For tangled magnetic fields however, a tangential 
component should be found. The criterion to suppress the KH 
instability is similar to before. It is suppressed when 



2n fi Picm Vj 2 CM \ 5fiG J 

( picm fV vicm r 2 do) 

\10- 27 gcm- 3 / \ 1000 km s- 1 / V ' 

implying that the KH instability occurs in weaker winds, even 
if the RT instability does not take effect, apart from very weak 
winds such as those in the cluster outskirts. 



3.2.1. Self-gravity 

We consider the extent to which self-gravity of stripped gas 
clouds can prevent RT instability. This case will apply at the 
cloud's surface, if the gravitational acceleration exceeds 
(Eq. HJ. If we assume a spherical gas cloud of density, p c i, and 
radius, R c \, after completing some simple arithmetic, we find that 
the RT instability is suppressed if 

Ed n c i R c \ 



10 21 



Picm 



JO" 27 gem- 3 



cm" 3 kpc 



VlCM 



1000 km s- 



(8) 



This implies that only clouds of high column densities can pre- 
vent an RT instability. We remark, however, that these clouds 
may yet be affected and dissolved, by the KH instability. 

3.2.2. Magnetic fields 

A second mechanism to stabilise stripped-off clouds against RT 
instability are tangential magnetic fields. We assume again a 
spherical gas cloud of radius, R c \, and substitute the expression 
for S i, from Eq.[2]into Eq.[5]to obtain 



A r = 



B 2 cos 2 9 f fid 

i 2 
po jCw Picm v icm 



= 1 .06 ■ /? c i cos 2 ^^) 



Picm 



10~ 27 g cm 



vicm 



1000 km s" 



■rm 



(9) 



This result is independent of the density of the cloud, and com- 
pares the energy density of the ICM wind with the energy density 
of the magnetic field. The magnetic field in spiral galaxies con- 
sists of an ordered large-scale field of approximately lfiG, and 
a small-sc ale, tangled co mponent with a field strength of about 
5nG (see iLequeuxl 120051 and references therein). We therefore 
assume that because of the small-scale field, a magnetic field B- 
component tangential to the cloud's surface with field strength 
of the order of a few pG will exist. The critical value of A c is 
then of the same order as the radius of the cloud. It is therefore 
likely that the stripped gas will suffer RT instability only in case 
of strong ICM winds. 

3.2.3. Kelvin-Helmholtz-instability 

In addition to the RT instability, the stripped gas can suffer KH 
instability. This instability cannot be prevented by a gravitational 



4. Discussion 

In this paper, we have derived two main results. We have shown 
that the inner parts of ram pressure stripped disk galaxies, which 
are resistant to RPS, are in addition stable against RT instability 
in the case of face-on stripping. The same should be true at mild 
inclinations, as neither the gravitational acceleration in the wind 
direction, nor the projected gas surface-density changes signifi- 
cantly with inclination angle. In addition, the clumps of stripped- 
off gas are exposed to RT instability and should be destroyed. 

Given that the RT instability is suppressed for the inner gas 
disk that survives ram pressure stripping, we conclude that the 
structures observed at the leading edge of NGC 4402 are not 
caused by the RT instability. These structures could be caused 
by different processes, one possibility being the inhomogeneous 
ISM as discussed by Crowl et al. (2005). Quilis et al. (2000) first 
studied the effect of an inhomogeneous ISM on ram-pressure 
stripping simulations, by creating holes inside the gas disk. The 
holes produced more pronounced Kelvin-Helmholtz instabili- 
ties, and thus stronger mass loss. The observations of NGC 4402 
do not concern under-dense but over-dense regions in the ISM: 
dense gas clumps appear to be left behind after the diffuse gas is 
stripped, but the dense clouds also appear to be ablated. To fol- 
low such structures in simulations requires extremely high reso- 
lution and computational effort. Nonetheless, this is an important 
task for future simulations. 

To what extent could stripped gas be observed? A com- 
plete investigation of this question would have to include pro- 
cesses such as thermal conduction and evaporation, cooling 
and star formation, as well as the effects of RT and KH in- 
stabilities. The above estimates imply that the stripped gas is 
likely affected by RT instability, and possibly also KH insta- 
bility. Nevertheless, more massive stripped-off gas packages 
can be gravitationally bound and stable against RT instability. 
Qualitatively, this can explain why only a few massive clouds are 
visible (in H ff ) for g alaxies observed being affec ted by RPS (such 
as e.g. N GC 4522: kennev & KoopmannlT999t and NGC 4569: 



Fig. 4 in lTschoeke et al.ll2001l) while in smaller clumps the de- 
structive RT instability dominates and disperses the gas into the 
ICM. The fact that the stable, stripped-off clouds emit H a in 
spite of no stars having formed, implies that they are heated 
by external sources. As the most plausible candidate, heat con- 
duction should be suggested because, in contrast to general as- 
sumptions, self-gravitating clouds are not de stroyed by evapo- 
ration but accumulate gas by condensation dVieser & Henslerl 
1 20071). and ar e, furt hermore, stabilised against KH instability. 
Vollm er et al. discussed the influence of thermal conduc- 

tion and effect of reduced heating on stripped gas clouds. They 
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argued that the clouds no longer absorb the stellar far-ultraviolet 
radiation field. They conclude that the stripped gas could be ei- 
ther ionised and hot, or in the form of cool and even molecular 
clouds, enclosed in neutral atomic shells. More detailed studies 
are required to study the combined influence of these effects, as 
well as the influence of star formation. 
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